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Abstract 

The (cr, uS) model in mean-field approximation where the meson fields are treated 
classically, describes much of observed nuclear structure and has been employed 
to describe the nuclear equation of state up to the quark-gluon phase transition. 
The acceleration of the meson sources, for example, in relativistic heavy-ion colli- 
sions, should result in bremsstrahlung-like radiation of the meson fields. The many 
mesons emitted serve to justify the use of classical meson fields. The slowing of the 
nuclei during the collision is modeled here as a smooth transition from initial to 
final velocity. Under ultra-relativistic conditions, vector radiation dominates. The 
angular distribution of energy flux shows a characteristic shape. It appears that if 
the vector meson field couples to the conserved baryon current, independent of the 
baryonic degrees of freedom, this mechanism will contribute to the radiation seen 
in relativistic heavy-ion collisions. The possible influence of the quark-gluon plasma 
is also considered. 

Key words: relativistic heavy-ion collisions; vector meson production; 
bremsstrahlung; quantum hadrodynamics; relativistic mean-field theory 
PACS: 25.75.-q, 24.10.Jv, 13.60.Le 



1 Introduction 



The (cr, a;) model of the nucleus [1,2] is a relativistic quantum field theory 
which describes the nuclear interaction using three fields. They are neutral 
scalar meson, neutral vector meson and baryon fields. The scalar meson field 
couples to the scalar density, ^t/), while the vector meson field couples to the 
conserved baryon current, = ip'-f^ip. In the relativistic mean-field approxi- 
mation (RMFT) where the sources are large, the meson fields can be replaced 
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by classical fields and the linearized Dirac equation solved exactly. This model 
has had several successes, properly describing: bulk properties of nuclear mat- 
ter, ground state properties of nuclei, the excitation spectrum of nuclei, low 
energy nucleon-nucleus scattering observables, as well as collective motion. All 
except the first of the preceding list incorporate spatial dependence for the me- 
son fields and their sources. The last item, collective motion, incorporates a 
slow time dependence in addition to the spatial dependence. 

It is here proposed to use this model to calculate vector meson production 
during relativistic heavy-ion collisions like those to be seen at RHIC. All of 
the above successes of the (a, u) model have been for low and intermediate 
energy phenomena. Therefore, some justification is necessary for use of this 
model in a radically different physical situation. First, using ground state wave- 
functions calculated from this model and empirical nucleon-nucleon scattering 
amplitudes one can obtain excellent agreement with experimental nucleon- 
nucleus scattering observables in the relativistic impulse approximation up to 
energies of the order 1 GeV [3]; second, the longitudinal response for 56 Fe(e,e') 
at momentum transfer of |<f| = 0.55 GeV has been calculated by Frank [4] in 
this model and agrees well with data; in addition, the (a, u) model taken in 
RMFT has been employed to describe the nuclear equation of state up to the 
quark-gluon phase transition [1,5]. Furthermore, the presence of many quanta 
of the meson fields here serves to validate the classical approximation. 

It is not clear what the appropriate degrees of freedom actually are for the 
RHIC-like collisions.^ Although the (a, to) model deals with hadrons, these 
results will turn out to be more general. We consider a process whereby two 
conserved baryon currents pass through each other, are slowed, and thus ra- 
diate energy via bremsstrahlung in the form of vector mesons. In the end, the 
only requirements for this calculation are an empirical knowledge of the initial 
and final rapidities of the baryon currents, conservation of these currents, and 
a vector meson that couples to the conserved baryon currents. In the final 
analysis, it is irrelevant whether the baryon current is carried by nucleons or 
quarks during the collision, as long as the vector meson couples to the baryon 
current. 

It should be noted that in the work of Mishustin et. al [6] , where the (cr, to) 
model is employed to calculate baryon-anti-baryon production via virtual vec- 
tor meson bremsstrahlung under RHIC conditions, it is suggested that real uj- 
mesons should be produced via bremsstrahlung and that this process should 
be relatively soft.^J 

2 In the work of Mishustin et. al [6] it is stated, "One cannot even say what degrees 
of freedom, hadrons or quarks and gluons, are more suitable for describing these 
collisions." 

3 In Ref. [6] it is further observed about vector meson bremsstrahlung that "... it 
is clear that the same mechanism can produce also mesons. For instance, the real 
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In this work only central collisions of identical heavy nuclei are considered. 
Therefore, the lab and equal- velocity frames are equivalent. In Ref. [7] Hanson 
shows evidence that for central high-energy collisions (rapidity larger than 3), 
the nuclei are expected to be slowed by the collision, but not completely 
stopped. This slowing occurs over the short time of the collision. This gives 
rise to a large deceleration. In this model, the decelerating sources radiate the 
classical meson fields. This is the same process as classical electromagnetic 
bremsstrahlung for an accelerated charge. 

Due to Lorentz contraction, the scalar meson radiation is greatly suppressed 
relative to the vector meson radiation. Most of the vector meson radiation 
is in the form of high energy mesons whose mass is therefore neglected. An 
attenuation factor is included to account for the strong interaction between 
the vector mesons and the baryons. The attenuation factor is constructed in a 
Lorentz invariant manner to maintain the covariance of the theory. This model 
contains three adjustable parameters: the deceleration time, the rapidity loss 
and the total cross section for the vector meson-nucleon interaction. The an- 
gular distribution of energy flux is characteristic, varying only in magnitude 
for variations in the cross section and deceleration time. Variations in rapidity 
loss affect both the magnitude and width of the angular distribution; however, 
the general shape is robust for parameter variations. 

This work indicates that the bremsstrahlung of vector mesons could contribute 
significantly to the total radiated energy during relativistic heavy-ion colli- 
sions. The model used here is very simple; it depends only on the vector mesons 
coupling to a conserved baryon current. The fact that this model predicts an 
appreciable amount of radiation with a characteristic angular distribution that 
is robust against parameter variation suggests that it is deserving of further 
investigation. 

This radiation of mesons during relativistic heavy-ion collisions has been con- 
sidered before. Weber et. al [8] examined the dynamics of relativistic heavy- 
ion collisions using the (a, uo) model in mean field theory. They modeled the 
collisions within the relativistic Boltzmann-Uehling-Uhlenbeck model using 
full solutions to the meson- field equations. At the energies they studied (1- 
20 GeV/A), they found meson radiation to be negligible (< 3 MeV/nucleon). 
However, the amount of energy radiated increases as the acceleration squared, 
so their work does not rule out appreciable radiation at higher energies. The 
production of pions and photons via bremsstrahlung was investigated in Refs. 
[9,10]. Ivanov [11] calculated a>meson radiation due to filamentation instabil- 

w-meson can be generated in the bremsstrahlung process when the four-momenta 
of quanta satisfy the mass shell constraint p 2 = m^. ... These channels are charac- 
terized by lower threshold and, therefore, by smaller momentum transfers. Since the 
corresponding coupling constants are also large, one can expect high multiplicities 
of mesons coherently produced in ultra-relativistic nuclear collisions." 
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ity. The behavior of the meson fields during relativistic heavy-ion collisions 
was explored in Refs. [12,13]; however, these are at much lower bombarding 
energies and did not directly look at meson radiation. 

It is expected that at the energies available at RHIC these central collisions 
will produce a quark-gluon plasma. The question then arises, what effect will 
the creation of the plasma have on this model that uses baryons and mesons? 
Indeed, would any sign of this hadronic process survive under these condi- 
tions? It has been suggested that as individual quarks interact during the 
collision, color strings or flux tubes will form [14-16]. These flux tubes are 
stretched out behind the nuclei as they pass through each other (see, for ex- 
ample, figure 4 of Ref. [14]). The breaking of these flux tubes then provides 
an important contribution to the formation of the quark-gluon plasma. The 
meson production described in the present work comes exclusively from the 
forward going baryons that constitute the initial colliding nuclei. Even if the 
mesons must pass through part of the baryon-rich quark-gluon plasma, there 
is theoretical evidence that meson-like modes can propagate through a quark- 
gluon plasma [17]. Such a physical situation may alter the attenuation but 
could still leave a detectable sign of this mechanism. 

This paper attempts to test the limits of the (a, to) model in RMFT for a 
strong time-dependence. Section 2 introduces the framework for the (a, uj) 
model and classical bremsstrahlung. Sections 3 to 5 discuss the model for the 
baryon current, the degree of incoherence in the radiation and the attenuation 
of the vector meson radiation. The energy spectrum is investigated in section 6. 
Sections 7 and 8 discuss the results and conclusions. 



2 Formalism 

The basic (cr, to) model is defined by the lagrangian density,^] 



where is the vector field, <p is the scalar field and ip is the baryon field. The 
masses of the nucleon, vector meson and scalar meson are Mat, and m a 

4 Throughout this paper natural units are used, c = % = 1. In these units 0.197 
GeV = 1 fm . The metric is that of Bjorken and Drell and Ref. [2]. 









(1) 
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respectively. The field tensor is defined as 

w dvu 2 
dx^ dx u 

This lagrangian neglects non-linear self-couplings of the scalar field. In the nu- 
clear ground state, the vector and scalar fields are of comparable strength [1]. 
Under a Lorentz boost the scalar density is invariant; however, the vector 
density is enhanced by a factor of 7, defined as 

7 = (l-/5 2 )- 1 / 2 , ( 3 ) 



with (3 the velocity. For the dynamics of the collisions investigated in this pa- 
per, this factor varies between 10 and 100. The radiated energy is proportional 
to the field squared; therefore, the vector meson radiation is at least a factor 
of 100 greater than the scalar meson radiation. Consequently, the scalar field 
can be neglected. 

The energy-momentum tensor is, 

fir da 1 8V a 



This form of the energy-momentum tensor works well for calculating the total 
energy and momentum of the system. However, it contains total divergences 
that will give incorrect results for the energy flux. Therefore, the symmetric 
energy-momentum tensor must be used [18], 

9^ = ^ v F af5 F a p + F^F/. (5) 



The power radiatedP] is 

dP(t') 



an 



R 2 [S(t)-n] ret (6) 



where P(t') is the power radiated at time t' and <%(= 0° l ) is the Poynting 
vector. The observer is located a distance R from the source of the radiation 
along the unit vector n. The notation, [] ret , means evaluate at t' — t — R(t'). 

The equation of motion for the vector meson field is 

8 V F^ + rnlV* = g v B». (7) 



The development of the radiated energy follows closely that in Jackson [19]. 
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Since the baryon current is conserved, this reduces to [d 2 + m^V 1 = g v B^. 
This is the inhomogeneous Klein-Gordon equation for a massive particle which 
has the solution, 



In the above equation, x is the four-position (t,x), k is the four-momentum 
(e, k) and B^(k) is the Fourier transform of the current. A convergence factor, 
i?7£, is introduced to eliminate the singularity at k 2 = m 2 . The choice of — vqe 
gives the retarded solution. The solution above can be converted to an integral 
over d 4 x' and gives V tx {x) = J d 4 x'Q r (x — x')B^(x') where 



(27r) 4 J m 2 — k 2 — ir/e 



(x-x') = -—- [ d A k 
y ! 2tt 4 J 



8(t - Tl) m. 



4vr^ 4ttVt 2 - K 2 



Jiim^Vr 2 -n 2 )&(T-n), (9) 



where r = t — t' and 1Z — \x — x'\. J\ is the first order Bessel function and 
Q(y) is the Heavyside step function defined as 0(y) = (1 + \y\/y)/2. 

If the baryon current is replaced by the electromagnetic current, the coupling 
constant changed to the fine structure constant, g 2 /^ — > a, and the mass, 
m w , goes to zero, this reduces to electromagnetic bremsstrahlung. Note that as 
rriu — > the second term in equation (9) vanishes. It is argued in section 6 that 
neglecting the mass of the vector meson does not alter the results appreciably, 
and we will take = for the rest of this work. 

The above solution for and equation (5) gives a distribution of energy as 
a function of outgoing energy, e, and angle of, (compare equation (14.70) of 
Ref. [19]) 



d 2 E e 2 g 2 



de dtt 4tt 2 



OO 

J dt J d 3 xn x [n x B(x,t)\ e 



ie(t—n-x) 



(10) 



Again, n is the unit vector from the point of radiation to the observer. 

The radiation is strongly peaked in the direction of the baryon cluster's mo- 
tion. Thus, we concentrate now on just one of the incident nuclei. Assuming 
the current is that of a collection of point baryons traveling with identical 
velocity, 

B(x,t)=J25(x-x J (t))(3(t) (11) 

j 
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we may do the spatial integration. The baryon's position is given as Xj(t) = 

— * — * 

X(t) + fj(t) where X(t) is the position of the center-of-mass and fj(t) is the 
baryon's position relative to the center-of-mass. The relative position of the 
baryon will vary with time, for example, due to changing Lorentz contraction 
and the empirically seen spreading of final rapidity; however, in this simple 
model this motion will be neglected relative to the much more important 
center-of-mass motion. The baryons will this be considered as "frozen" relative 
to the baryon clusters as they pass through each other and decelerate. The 
baryon's position thus becomes, 



x 



t 

;{t) = J P(t') dt' + fj(-00) = X(t) + fj. 



(12) 



Separating out the dependence on j we may write, 
d 2 E 



e 2 gl 



de dil 4tt 2 



dtn x [fix (3(t)\ e 



ie(t-n-X(t)) 



(13) 



Taking the sum outside the square gives a coherence factor 

A A 



j=i k=i 



-ien-(rj-r k ) 



(14) 



One must now take a statistical average over the positions of the baryons in 
the incident nucleus; this average is weighted by the square of the ground state 
wave-function. If the exponent is large, as argued in section 4, the off-diagonal 
elements average to zero and the coherence factor is approximately A, the 
number of baryons. 

We now wish to find the angular distribution of energy flux and the total en- 
ergy radiated. Equation (13) must therefore be integrated over all frequencies. 
Equation (13) may be rewritten as 



d 2 E 



9l 



de dVl 4tt 2 



A 



dt 



fix [fix (3(t)] c i £( t-n-x(t)) 



(15) 



where we have assumed the velocity and acceleration are collinear to write 



dt 



n x [n x P(t)] 
1-n- (3(t) 



fix [fix $(t)} 



(16) 
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Changing variables to r = t — ft ■ X(t) and expanding the square we may now 
do the e integration. The integration should be over positive energies; however, 
the integrand is even in e so we take 1/2 the integral over all energies. The 
angular distribution of energy flux is 



dE 



9 i A S 



df 



n x [n x ${t)} 



[l-n-Pit)f 



(17) 



The angular integration can also be done and this gives 



E = ^A I dt 7 6 |/5(t)| 2 . (18) 



Note that the last two results are just the single particle result (see equations 
(14.38) and (14.43) of Ref. [19]) times the coherence factor, V « A. 



3 Baryon Current Models 



The only quantity required to evaluate equation (10) is the baryon current, 
Bix,t). In the present covariant approach, this is related to the particles' 
velocity, /3, in the lab frame by 

B,ix,t)= £ P { i\x)U^it) (19) 

A=l,2 



where 1 and 2 are the two nuclei, pb is the rest-frame baryon density and = 
(7, 7/3) is the four-velocity of the nucleus. The baryon current is conserved; 

° = 0, (20) 



dx 



and transforms as a four-vector, A^B^Ax) = X)a=i,2 Pb \Ax)A^lJ( X) iAx). 

The baryon density should be determined from the ground-state nuclear wave- 
functions, computed in the Dirac-Hartree approximation, summed over the 
occupied orbitals. This procedure is explained in detail in Ref. [20] and im- 
plemented in Ref. [21]; however, it is argued in sections 2 and 4 of this paper 
that the nucleons radiate incoherently and the energy flux is insensitive to the 
exact baryon distribution. Therefore, one can approximate the true baryon 
density by a constant density inside the nuclear radius, ps(r) = p o 0(i?jv — r). 



8 



The motion of the baryons will be modeled as having a smooth transition from 
their initial speed to their final speed. A Fermi-type parameterization of the 
speed is used [6]: 

m = P f + ^^, (2i) 

where r is the stopping parameter. Call the unit vector along the beam axis z. 
Therefore, the first nucleus is traveling with velocity j3z and the second with 
—/3z. The initial and final speeds are fixed from experimental results and the 
stopping parameter is allowed to vary in this model. 

The beam energy per nucleon is E heSim /A. The initial rapidity is related to the 
beam energy by M N coshyi = E heiim /A, and the initial speed by fa = tanhyj. 
The energy lost during the collision is characterized by the rapidity loss, Sy = 
yi — yf. Mishustin et. al [6] claim a rapidity loss of by = 2.4 ± 0.2 for central 
Au+Au collisions at RHIC energies of E^, eam /A = 100 GeV. The same value 
is used for this study of Pb+Pb at the same energy. The stopping parameter 
is allowed to vary between 5 and 20 fm.0 



4 Incoherence 

The initial speed of the nucleons for E hea _ m /A = 100 GeV is fa ~ 0.99996. The 
velocity of the nucleons relative to the center-of-mass is small compared to 
the velocity of the center-of-mass. This means that the collision time between 
nuclei will be much shorter than the interaction time between nucleons in the 
same nucleus. The nuclei then can be treated as static collections of nucleons 
with uniform longitudinal motion. This is similar to the frozen approximation 
used in the parton model. 

The baryon current can be represented as a collection of point particles: 

A A 

B(x, t) = pW(t) £ 5[x - X x {t) - r 3 ] + p®(t) £ 5[x - X 2 (t) - fj] (22) 

3=1 j=l 

where the first (second) term represents the first (second) nucleus, Aj(i) is the 
position of the center of mass of the i th nucleus and fj is the position of the j th 
baryon relative to the center of mass at t = — oo. Since the lab frame coincides 
with the equal-velocity frame, X\ = = X. The position is the integral of 

6 Investigations up to 50 fm still maintain the characteristic shape of the angular 
energy flux distribution. 
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the velocity with the boundary condition that the velocity asymptote to the 
free particle solution at t — ±00. The exact form used is 



X(t) = J pit') dt' z = frtz + i(3 f - Pi)r ln(l + e t/r )z. (23) 
Using this form of the current in equation (10) gives 



^ = I? V I dt I dt ' [m ***WW) sin ^ ei£{t ~ t,] ( 24 ) 

x sm[e(f>it)} sm[e<f>(t')] 

where sin^(i) = n x (n x z) at time t, 0(t) = n • and V is the coherence 
factor defined in equation (14). 

If in the coherence factor, \en ■ (fj — f k )\ <C 1 then the phases are all near 
zero. If, on the other hand, \en ■ (fj — fk)\ > 1 then the phases are large and 
nearly random. In the end, when a statistical average over all possible con- 
figurations is taken using the ground-state wave-functions, the first case leads 
to V ps A 2 corresponding to the nucleons radiating as one large baryon of 
"baryonic charge" A. For the latter case the coherence factor becomes V ^ A 
corresponding to each nucleon radiating independent of its neighbors. Since 
neither of these conditions hold absolutely, the coherence factor will be some 
power of the atomic weight, A a where 1 < a < 2. To neglect the vector meson 
mass, its outgoing energy must be on the order of several GeV. The inter- 
particle distance is on the order of 1 fm. The latter condition holds and the 
nucleons should radiate incoherently, a — 1. 



5 Attenuation 



Throughout this work the vector meson field has been treated analogously to 
the electromagnetic field with just a change in coupling strength. However, 
there is one aspect where there is a significant difference. The vector mesons 
interact strongly with the baryons; therefore, there will be an attenuation 
of vector mesons as they pass through either the nucleus from which they 
are radiated or the other nucleus. The attenuated energy flux can appear 
as radiation through another channel or could be used to heat the nuclear 
material. Since most of the radiation will occur early during the collision 
[equation (17) is dominated by the denominator which increases during the 
collision] and is directed in the forward direction, this attenuation is non- 
negligible. This reduction is accounted for by including a multiplicative factor, 
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A 2 {t), in the radiated power, 



dP(tf) 



R 2 [A 2 (t)S(t) -n) 



(25) 



dfi 



Forms of equations (17) and (18) involving the attenuation factor follow the 
development in section 2. 

This attenuation factor[[] is taken as, 



where I and x are the distances traveled through the radiating nucleus and 
the second nucleus respectively, a is the vector meson-nucleon total cross- 
section, and pb is the baryon density. Equation (26) can be evaluated in the 
rest frame of the radiating nucleus making I time- independent. To maintain 
the covariance of the model the attenuation factor is evaluated along the 
axis of motion and taken to be independent of outgoing angle. Since most 
of the radiation is in the far forward direction, this results in only a small 
overestimation of the attenuation. 

To evaluate equation (26), consider the radiation coming from the nucleus 
traveling in the positive z direction (the argument is identical for the other 
nucleus). Move the origin for the integration to the center of mass in nucleus' 
rest frame. The distance traveled inside this nucleus depends only on where 
the radiation originates (remember, the attenuation is assumed independent 
of outgoing angle). Using cylindrical coordinates (p, (ft, z) to take advantage of 
the azimuthal symmetry, we find 



where Rn is the nuclear radius. Since the radiation must come from inside the 
nucleus the following condition always holds: z 2 + p 2 < R 2 N . The form of the 
equation for the distance through the target nucleus is more involved since 
the nucleus is moving. Defining the distance between the two centers of mass 



It is assumed here that the vector meson will interact with the baryons immedi- 
ately after it is created. If the nucleus demonstrates color transparency, then this 
attenuation factor will be reduced. 



A(t) 




(26) 




(27) 



11 



as d — —2j'(t)X(t), the distance traveled through the moving nucleus is 



x(t) 



-( z -d) + ^R%-p? 



d — Rn /l' > z 

d- R N /i < z <d + R N /i (28) 
otherwise 



where 7' is the Lorentz factor for the target moving in the rest frame of 
the radiating nucleus. It is related to the Lorentz factor in the lab frame 
[equation (3)] by 



7 / (t)=7 2 (t)(l + /? 2 



1 + P 2 
l-(3 2 



(29) 



6 Energy Spectrum 



Throughout this paper the vector meson mass has been assumed small com- 
pared to the particles' outgoing energy, and hence neglected. It is the purpose 
of this section to justify this assumption. Equation (24) gives the energy dis- 
tribution of the radiated mesons. If the observer is assumed far away, the 
angle to the observer, ip(t), will change little over the time of interest. The 
sin^ factor may be taken outside the time integral as a constant. We define 
the meson energy spectrum by keeping only the energy and time dependent 
factors in equation (24), 

2 

(30) 



X\e) 



dt f3(t) e iet sm[en-X(t)} 



Note that X 2 ^) is dimensionless. Again assuming the angle to the observer 
changes little and the velocity remains near one, the integral can be evaluated 
analytically giving: 



l 2 {e) = {erf 



B[ier(l — f3i cos ip), ier(f3f cosip — 1)] 



+ B[ier(l + $ cos^), —ier(f3f cosip + 1) 



(31) 



where B[x,y] is the beta function which is related to the gamma or factorial 
function by 



B[x,y] = 



T(x + y)' 



(32) 
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1 ■ 1 ■ 1 ■ 1 ■ 1 

50 100 150 200 

e [GeV] 

Fig. 1. The dimensionless energy spectrum, X 2 (e), as denned in equation 30 as a 
function of outgoing vector meson energy. 

As can be seen in figure 1 the peak occurs at e = 0, and decreases slowly as 
e increases. This slow decrease means that most of the radiation is above the 
vector meson mass and that the massless assumption is acceptable. Cutting the 
spectrum off at the meson mass e = = .783 GeV, reduces the integrated 
spectrum by less than 1%. 



7 Results 

Throughout this work, only central collisions of identical nuclei at a bombard- 
ing energy achievable at RHIC of E^ eam /A = 100 GeV [22] have been consid- 
ered. A short discussion is still needed concerning the value of the required 
coupling constant, g v . There is no solid evidence to direct a choice of the cou- 
pling constant for the type of reaction that this work describes. A value that 
reproduces static nuclear properties has been chosen [1]; <?^/47r = 10.8. Using 
a value of g^j^K between eight and nine, Gross, Van Orden and Holinde [23] 
were able to reproduce free nucleon-nucleon phase shifts to several hundred 
MeV. Since the {a,u) model is not an asymptotically-free theory, the cou- 
pling constant may even increase at higher energy and momentum transfer; 
this would lead to even more radiation than is shown below. The paper by 
Mishustin et. al [6], uses a coupling constant of g^/kit = 15.1, stronger than 
what is used here. As long as the coupling constant is of at least order 1, this 
model predicts appreciable amounts of energy radiated through this channel. 
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Fig. 2. Angular distribution of energy 
flux as a function of lab angle for out- 
going vector mesons. The rapidity loss 
is fixed at 5y = 2.4 and the vector me- 
son-nucleon total cross-section is fixed 
at a = 30 mb. It is shown for r = 5.0 
fm (solid line), r = 10 fm (dashed line) 
and r = 20 fm (dotted line). The A 
and g%/4ir factors have been divided 
out. 



Fig. 3. Same as figure 2 except stop- 
ping parameter is fixed at r = 10 fm 
and the cross-section takes the values 
a = 20 mb (solid line), a = 30 mb 
(dashed line) and a = 40 mb (dotted 
line). 



This is a crude model that has incorporated several simplifying assumptions. 
It would be inconsistent to think of the predicted numbers here as more than 
a rough guide to the order of magnitude of what may actually be seen. To 
remove any ambiguity in the coupling constant, the angular distributions are 
shown with the coupling constant divided out. 

Figures 2, 3 and 4 show the angular distribution of radiated energy in the 
form of vector mesons. Following the analysis of section 4, the energy flux is 
calculated as the incoherent sum of 2A nucleons radiating as point particles, 

dVl 4vr £ Mjl -/?(£) cos #) 5 (1 + /3(t) cos 9) 5 J ' 1 J 

The factor (g%/4ir)A is divided out in the figures to remove any ambiguity.^ 

8 The dimensions of the remaining expression are converted by 1 fm" 1 = 0.197 
GeV. 
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e [deg] 

Fig. 4. Same as figure 2 except stopping parameter is fixed at r = 10 fm and the 
rapidity loss takes the values Sy = 2.2 (solid line), Sy = 2.4 (dashed line) and 
5y = 2.6 (dotted line). 

In each of the figures, two of the parameters are held constant while the third 
is varied over reasonable values. The peak of energy flux occurs around 2.5° 
off the beam-axis. It is seen that variations in the stopping time, r, and the 
cross-section, a, only affect the magnitude and not the characteristic angular 
distribution shape. Changes in the rapidity loss, Sy, affect both the magnitude 
and shape of the distribution with larger rapidity losses in this range causing 
less radiation and causing the peak to be moved further off axis and broadened. 
It is interesting to note that the smallest rapidity loss used in this study, 
Sy = 2.2, results in the most energy being radiated. This can be explained by 
considering the circumstances under which most of the energy is radiated. At 
relativistic speeds, the energy radiated is dominated by the factor of (1— n-/3)~ 5 
in equation (17). A smaller rapidity loss results in the nucleus spending more 
time traveling faster, resulting in more radiation. Of course this only holds over 
a small region. At some point the smallness of the acceleration overcomes the 
smallness of the denominator and at zero rapidity lost, there is no radiation. 

For the intermediate values of the parameters (r = 10 fm, a = 30 mb and 
Sy = 2.4) and A = 208 (Pb beams) the total energy radiated, E rad , is 353 
GeV. The total energy available to be radiated in this collision is 3.78 x 10 4 
GeV. Reducing the stopping parameter to r = 5 fm gives _E ra d = 730 GeV. 
A rapidity loss of Sy = 2.2 and a stopping parameter of r = 10 fm gives 
-^rad = 516 GeV out of a total available energy of 3.70 x 10 4 GeV. On the 
order of 1% of the total energy loss is through this channel. 
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8 Conclusions 



In this paper a model for treating the bremsstrahlung radiation of neutral 
vector mesons coupled to a baryon current during central relativistic heavy- 
ion collisions has been developed. This model treats the nuclei as clusters of 
baryons frozen in relative position over the time of the collision. The clusters' 
velocities are modeled as changing smoothly with time. Lorentz contraction 
greatly increases the baryon density relative to the scalar density; hence, the 
scalar meson radiation becomes negligible at high energies. Most of the energy 
is radiated as highly energetic vector mesons, allowing for their mass to be 
neglected. The vector mesons interact strongly with the baryons and therefore, 
an attenuation factor is included. To maintain the covariance of the model, 
this attenuation factor is assumed independent of angle. 

The modeling of the flow of the baryon current is the only freedom in the 
model. By using a smooth connection between the initial rapidity of the nuclei 
with the experimentally measurable final rapidity, the baryon current depends 
on just two parameters, the rapidity loss (obtained from the initial and final 
rapidities) and a stopping parameter which is allowed to vary freely. Although 
this model is used well outside its tested domain, in the end only the vector 
field coupling to the conserved baryon current is seen. The model predicts a 
characteristic shape for the angular distribution of radiation that is robust 
against parameter variations. It appears this mechanism may contribute to 
meson production in the next generation of relativistic heavy-ion collisions 
like those at RHIC. This possibility would seem to be worth exploring exper- 
iment allyQ 
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If this process were to be observed, it could provide a diagnostic for the baryon 
flow during such collisions. For example, only straight line motion of the baryons 
has been considered. The radiation is focused sharply around the instantaneous 
motion of the baryon current so deviations from the angular distribution shown 
here could be used to track the intermediate flow of baryons during the collision. 
See also Ref. [11]. 
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